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Abstract. We have studied the Globular Cluster System of the merger galaxy NCC 1316 in Fornax, using CCD 
BVI photometry. A clear bimodality is not detected from the broadband colours. However, dividing the sample 
into red (presumably metal-rich) and blue (metal-poor) subpopulations a.t B — I — 1.75, we find that they follow 
strikingly different angular distributions. The red clusters show a strong correlation with the galaxy elongation, 
but the blue ones are circularly distributed. No systematic difference is seen in their radial profile and both are 
equally concentrated. We derive an astonishingly low Specific Freq uency for NGC 1316 of only Sn = 0.9, which 
confirms with a larger field a previous finding by Grillmair et al. (1999). Assuming a "normal" Sn of ~ 4 for 
early-type galaxies, we use stellar population synthesis models to estimate in 2 Gyr the age of this galaxy, if 
an intermediate-age population were to explain the low Sn we observ e. This value agrees with the luminosity- 
weighted mean age of NGC 1316 derived by Kuntschner & Davies (1998) and Mackie & Fabbiano ( 1998| ). By fitting 
functions to the Globular Cluster Luminosity Function (GCLF), we derived the following turnover magnitudes: 
B = 24.69 ± 0.15, V = 23.87 ± 0.20 and I = 22.72 ± 0.14. They confirm that NGC 1316, in spite of its outlying 
location, is at the same distance as the core of the Fornax cluster. 
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1. Introduction 

The analysis of globular cluster systems (GCSs) in ellipti- 
cal galaxies can have different motivations. One of them is 
to investigate the variety of GCS morphologies in relation 
to their host galaxy properties in order to gain insight into 
the formation of cluster systems (see Ashman & Zepf 1997 
and Harris 2000 for reviews). 

On the other hand, GCSs have been suc cessfu lly em- 
ployed as distance indicators (Whitmore et al. 1995 , Harris 
2000| ). This is particularly interesting if the host galaxy is 
simultaneously host for a type la supernova whose ab- 
solute luminosity can accordingly be determined, as has 
been the case for SN 1992A in NGC 1380 (Delia Valle et al. 



NGC 1316 (Fornax A), the brightest galaxy in the 
Fornax cluster, is also one of the closest and brightest ra- 
dio galaxies. Catalogued as a SO peculiar galaxy, it hosted 
two supernovae of type la (SN 1980D and SN 1981N). The 
core of the Fornax cluster is quite compact, which makes 
it a better extragalactic distance standard than the Virgo 
cluster (see Richtler et al. |1999D . However, NGC1316 is 
located relatively far out, 3? 7 away from the central giant 
elliptical NGC 1399. 



199C) 



199?) and SN 1994D in NGC 4526 (Drcnkhahn & Richtler 



However, it can happen that both aspects are equally 
interesting as with the target of the present contribution, 
NGC 1316. 



The early studies of Schweizer ( |1980| & |l98lD indi- 
cated that NGC 1316 is very probably a merger galaxy. 
Schweizer (1980) discovered the presence of several irreg- 
ular loops and tidal tails, as well as an inclined disk of 
ionized gas rotating much faster than the stellar spheroid. 
The age of the merger event was estimated to about 4-10*- 
2 • 10^ yrs. 

Studies of gal axy m ergers have shown (e.g. NGC 1275, 
Holtzmann et al. |l992|, NGC 7252, Whitmore et al. |l993, 
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NGC 4038/NGC 4039, Whitmore & Schweizer |1995D that 
bright massive star clusters can form in large numbers 
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during a merger event. But despite the strong evidence 
for a previous merger in NGC 1316, the only indication 
for cluster formation is that Grillmair et al. (1999, here- 



Table 1. Summary of the observations. 



after Gr99), could not see a turnover in the GCLF at the 
expected magnitude. They interpreted this finding as an 
indication for an enhanced formation of many less massive 
clusters, perhaps in connection with the merger. However, 
their HST study was restricted to the innermost region of 
the galaxy. 

In contrast to other galaxy mergers, where, presum- 
ably caused by the high star-formation rate (Larsen & 
I200C 



Richtler 1999 , |200C| ), the specific frequency of GCs in- 
creases, Gr99 found an unusually small total number of 
clusters relative to the luminosity of NGC 1316 {My ^ 
—22.8, adopting a distance modulus to Fornax of /i = 
31.35, Richtler et al. 



200C) 



There is also evidence from stellar population synthesis 
of integrated spectra that N GC 131 6 hosts you nger p opu- 
lations (Kuntschner & Davies 



1998 



Kuntschner |2000| ) and 
thus the question arises, whether the surprisingly low spe- 
cific frequency is caused by a high luminosity rather than 
by a small number of clusters. Goudfrooij et al. ( 2000 , 
hereafter GoOO) obtained spectra of 27 globular clusters 
and reported that the 3 brightest clusters have an age 
of about 3 Gyr, indicating a high star-formation activity 
3Gyr ago, presumably caused by the merger event. 

These findings and the hope for a good distance via the 
GCLF were the main motivation to do the present study 
of the GCS of NGC 1316 in a larger area than that of the 
HST study. As we will show, this galaxy resembles in many 
aspects the "old" merger galaxy NGC 5018, wh ose G CS 
has been investigated by Hilker & Kissler-Patig ( 1996 ). 

The paper is organized as follows: in Sect. |2 we dis- 
cuss the observations, the reduction procedure and the 
selection of cluster candidates. The photometric and mor- 
phological properties of the GCS are discussed in Sect. ||. 
Sect. ^ contains our findings concerning the Specific 
Frequency. We conclude this work with a general discus- 
sion in Sect, n 



2. Observations and Reduction 

The B, V and / images were obtained at the 3.6 m tele- 
scope at La Silla during the nights 29th and 30th of 
December, 1997 (dark moon), using the ESO Faint Object 
Spectrograph and Camera, EF0SC2. The field of view was 
5f6 X 5^6 with a scale of O'.'32/pixel. During the first night, 
short- and long-exposures in each filter were centred on 
the galaxy. In the second night, a background field located 
about 5' away from the centre of NGC 1316 was observed, 
overlapping by 1' the observations of the first night. In 
addition, several fields co ntaini ng standard stars from the 
Landolt catalog (Landolt 1992) were acquired in each fil- 

I 
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ter, as well as some short exposures of NGC 1316. Fig. 
shows the combined frames from both nights and Table 
summarises the observations. 





Date 


Filter 


Exposure (sec) 


seeing 


Dec 


29, 1997 


B 


4 X 600 


I'.'l 


Dec 


29, 1997 


V 


5 X 300 


I'.'O 


Dec 


29, 1997 


I 


6 X 300 


I'.'O 


Dec 


30, 1997 


B 


4 X 600 


1'.'3 


Dec 


30, 1997 


V 


3 X 600 


1'.'3 


Dec 


30, 1997 


I 


3 X 600 


1'.'2 




Fig. 1. The entire observed field in the B filter, combined 
from both nights. North is up, east is left. 



2.1. Reduction 

The frames were bias-subtracted and flat-fielded using 
standard IRAF procedures. Bad pixels were replaced by 
interpolation with neighbour pixels, using the task fixpix. 
The frames in each filter were then combined to remove 
the cosmic rays and to improve the S/N ratio. Small shifts 
were applied when necessary to avoid misalignment. 

To perform the photometry, a relatively flat sky back- 
ground must be provided. Therefore, the galaxy light 
needs to be subtracted. While in the case of early-type 
galaxies this is normally achi eved by fitting elliptical 
isophotes (see, e.g. Jedrzejewski 1987 ), NGC 1316 exhibits 
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Table 2. General parameters of the target galaxy, from de Vaucouleurs et al. (1991) and Poulam (1988). 
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name a (2000) 


5 (2000) I 


h 


type mv 


B-V Veloc. [km/s] 


NGC1316 03''22™41':6 


-37° 06' 10" 240? 16 


-56? 69 


(R')SAB(s)0 8.53 ±0.08 


0.86 1793 ± 12 



some substructures near its centre that can not be mod- 
eled in this way. We therefore applied a median filter, using 
the task median and a grid size of 35 x 35 pixels. 

After the subtraction of each combined image from 
its corresponding median frame, the object detection and 
PSF photometry was run using DAOPHOT under IRAF. 
The threshold level was set to 3cr, and typically 2500 to 
3000 objects were detected in each image. Many of these 
detections, particularly in the central parts of the galaxy, 
are spurious because of the remaining photon noise af- 
ter subtraction of the galaxy model (see Sect. 2^). Since 
globular clusters in the Fornax distance are not resolved, 
additional aperture corrections are not necessary. The ob- 
ject lists of the three frames were matched to remove 
spurious detections and to leave only objects detected in 
all three filters. This reduced the sample to 911 objects. 
We used the (x, y) coordinates of this list as input for 
the program SExtractor (Berlin & Arnouts 1996), which 
performed the classification of point-like and galaxy-like 
objects by means of a "stellarity index", calculated by 
a trained neural network from ellipticity and sharpness. 
Although SExtractor offers the possibility of performing 
aperture photometry, we stress that we used it only to get 
a better star-galaxy separation than is achievable with 
DAOPHOT alone. 

2.2. Calibration of the photometry 

A large number of standard frames were acquired in each 
band during the second night (the first night was not pho- 
tometric). The colours of the globular clusters are well in- 
side the range spanned by the standard stars. Aperture 
photometry with several annuli was performed for ev- 
ery standard star to construct a growth-curve (see, e.g. 
Stetson 199C). Only stars with well defined growth curves 
were selected. From these curves we chose an aperture ra- 
dius of 30 pixels (9'.'6) for the photometry of the standard 
stars. We then fit the instrumental magnitudes {nij) to 
the equations: 



Table 3. Colour and extinction coefficients for om pho- 
tometry according to (nl). 



mb = Abo + i\/b + ^bi ■{B-V)+A^,2-X 
mv = Avo + A-/v + ^vi -{V - I)+A^2-X 
mi = Aio + Ml + A-a ■ {V - I) + A^ ■ X 



(1) 
(2) 
(3) 



where Aj^ and Mj are the zeropoints and the standard 
apparent magnitudes, respectively. Aji and Aj2 are the 
colour and extinction coefhcients. X is the airmass. The 
fitted parameters are listed in Table |[ 

To check our photometry, we compared the magn itude 
and colours for NGC 1316 quoted by Poulain (1988) with 
our values using the short-time exposures. The mean ab- 
solute differences in magnitude and colour at five aper- 



filter 


Aji 


Aj2 


rms of the fit 


B 


-0.029 ±0.019 


0.220 ± 0.038 


0.025 


V 


-0.049 ±0.010 


0.111 ±0.017 


0.014 


I 


0.023 ±0.015 


0.046 ±0.029 


0.027 



ture radii between 22'.'9 and 86'.'6 are: (Ay) = 0.013, 
(As-y) = 0.007 and (Av-/) = 0.006 mag, well below the 
rms of the fit (see Table ||). We then defined 5 local stan- 
dard stars in the field of NGC 1316 to set the photometry 
of both nights in a consistent way. 

2.3. Selection criteria 

Several criteria have been applied to select cluster can- 
didates, according to colours, magnitude, photometric er- 
rors, stellarity index and projected position around the 
galaxy. We assume that the clusters are similar to the 
Milky Way clusters. Adopting an absolute turnover mag- 
nitude (TOM) oi V = —7.60 for the galactic clusters 
(Drenkhahn & Richtler |1999|, Ferrarese et al. |200(]| ) and 
11 = 31.35 (Richtler et al.^200g), the TOM of NGC 1316 is 
expected to be ^ ~ 23.7 mag, and the brightest clusters 
about V ~ 20. 

Although no reddening corrections are normally ap- 
plied when looking towards Fornax (Burstein & Heiles 
1982), we cannot restrict the colours of the clusters in 



NGC 1316 to match exactly the galactic ones. One point 
is the smaller sample of galactic clusters. Besides, we must 
allow for significant photometric errors of faint cluster can- 
didates. Fig. H shows a colour-magnitude diagram for all 
objects detected simultaneously in B, V and / in both 
nights (before the selection) , together with the cut-off val- 
ues adopted as criterium for this colour. As can be seen, 
the majority of objects have colours around V ~ I ~ 1.0. 
Objects bluer than 0.5 mag are very probably foreground 
stars. A fraction of the data points redder than V—I = 1.6 
are background galaxies. 

We also tested in our images the robustness of the 
"stellarity index" computed by SExtractor. This index 
ranges from 0.0 (galaxy) to 1.0 (star) and varies for the 
same object by about 0.2 when classifying under different 
seeing conditions, except for the brightest objects, which 
are clearly classified. By visual inspection of the images, 
we are quite confident that bright galaxies are always given 
indices near to zero. However, this classification becomes 
progressively more difficult with fainter sources. 
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Fig. 2. A colour-magnitude diagram for all objects de- 
tected in B, V and / in both nights before the selection 
criteria (911 points) 

. The dashed lines indicates the selection criterium in the V — I 
colour (see Sect. E^.) 
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Fig. 3. The stellarity index computed by SExtractor as 
function of the V magnitude. The dashed lined indicates 
the cut value adopted in the selection criteria to reject 
background galaxies. 



Fig. H shows the "stellarity index" for the 911 objects 
in our sample (before the selection criteria) as a function 
of the magnitude. Two groups of objects having indices 
of ~ 0.0 and 1.0 can be seen, but it is apparent that 
faint clusters cannot be unambiguously distinguished from 
background galaxies due to the uncertainty of the stellar- 
ity index in the case of faint sources. Very similar results 
were obtained with our "artificial stars" (see Sect. |2.5D , 
where indices down to 0.2 were measured for faint objects 
that are constructed using the PSF model and, therefore, 
are expected to have a stellarity index of ~ 1.0. 

Guided by our experience with artificial stars, we de- 
fined the cut-off value for the stellarity index to be 0.35. 
Remaining galaxies will be statistically subtracted be- 
cause the same criteria are applied to the background field 
(see Sect. U). 



Finally, we rejected objects with photometric error 
larger than 0.15 mag. To summarise, the following criteria 
were applied to our sample of 911 sources detected in B, 
V and /: 

i) > 19.5 

ii) OA<B-V < 1.4 

iii) 0.3 < V"-/ < 1.8 

iv) stellarity index > 0.35 

v) error(y), error(i? — V), error(y — /) < 0.15 

375 objects met this set of criteria and are our globular 
cluster candidates. 



2.4. Background correction 

After the selection, there might still be some contamina- 
tion by foreground stars and background galaxies in our 
sample. To subtract them statistically, one needs to ob- 
serve a nearby field, where no clusters are expected, and 
to apply the same detection and selection criteria as with 
the galaxy frame. However, the analysis of our background 
field still shows a concentration towards the galaxy centre, 
which means that only part of the field can be considered 
as background. 

We used the radial profile of the GC surface densi- 



ties (see Sect. 3.3) to select all objects on the flat part of 
the profile, i.e., where the number of globular cluster can- 
didates per area unit exhibits no gradient. Fig. |ll| (top) 
shows that this occurs at r « 300", where r is the distance 
from the optical centre. Thus, all objects with galactocen- 
tric distances larger than 300" and matching the above 
criteria, constitute our background sample. 

We constructed a semi-empirical luminosity function 
of the background clusters with a technique described by 
Seeker & Harris ( 1993| ), where a Gaussian is set over each 
data point, centred at the corresponding observed magni- 
tude. The sum of all Gaussians gives a good representation 
of the background, without introducing an artificial undu- 
lation or loss of information due to the binning process. 
Fig. ^ shows the histogram of the background objects, us- 
ing a bin size of 0.4 mag, and the adopted semi-empirical 
function, which we use as backgrou nd i n the calculation 
of the luminosity function (see Sect. 3^) and the specific 
frequency. Admittedly, the numbers are small. The dip at 
V — 23 may be simply a result of bad statistics. On the 
other hand, as Table 6 shows, the background counts are 
small compared to the clusters counts. Therefore, errors of 
the order of the statistically expected uncertainty do not 
significantly influence our results and are accounted for in 
the uncertainties of the total counts in Table 6. 

2.5. Completeness correction 

To correct statistically for completeness, one needs to de- 
termine which fraction of objects are actually detected 
by the photometry routines. These 'artificial stars ex- 
periments' were performed using the task addstar in 
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Fig. 4. The semi-empirical luminosity function for the 
background field (dashed line). A histogram with a bin 
size of 0.4 mag is over-plotted for comparison. 

DAOPHOT. In one step, 100 stars were added in the sci- 
ence frames from 1/ = 20 to y = 25 in steps of 0.1 mag, 
distributed randomly to preserve the aspect of the image 
and not to introduce crowding as an additional param- 
eter, but with the same coordinates in the -B, V and / 
frames every loop. The colours of the stars were forced 
to be constant B — V = 0.7 and V — I — 1.0, close to 
the mean colours of the globular clusters in NGC 1316 
(see Sect. 3T). The object detection, photometry, classi- 
fication and selection criteria were applied in exactly the 
same way as for the globular cluster candidates. To get 
sufficiently good statistics, we repeated the whole proce- 
dure 10 times, using different random positions in each of 
them and averaging the results. In all, 50000 stars were 
added. 

However, the completeness is not only a function of 
the magnitude. Due to the remaining noise after galaxy 
subtraction, the probability of detecting an object near the 
centre of the galaxy is smaller than in the outer parts. We 
divided our sample of artificial stars into elliptical rings, 
in the sa me w ay as we did in deriving the radial profile 
(see Sect. 3^). The results of the completeness tests are 
summarised in Fig. ||. It can be seen that the 50% limit 
goes deeper with increasing galactocentric distance. 

3. Photometric and morphological properties 

3.1. Colour distribution 

In this section we discuss the colours of our GC candidates. 
No in terstellar reddening is assumed (Burstein & Heiles 



1982 ) and therefore, only internal reddening might affect 
the colours. 

The histograms in Fig. ^ show the colour distribution 
of the cluster candidates. A fit of a Gaussian function with 
free dispersions returned the values listed in Table ^ and 
is overplotted for comparison. 

As can be seen, the dispersion in the B — V and V — I 
histograms is completely explained by the photometric er- 
rors alone. This is not the case for the B — I histogram, 
due probably to the greater metallicity sensitivity of this 




Fig. 5. Completeness factors in four elliptical annuli. The 
probability of detection strongly decreases near the centre 
of NGC 1316. 

Table 4. Fit of a Gaussian function to the colour his- 
tograms of the GC candidates. 



colour 


centre 




B-V 


0.80 ±0.02 


0.13 ±0.04 


V-I 


0.95 ±0.01 


0.17 ±0.02 


B~ I 


1.77 ±0.02 


0.35 ±0.03 



colour. A Gaussian with a — 0.15 mag (the maximum al- 
lowed photometric error, see Sect. 2.3) is overplotted and 
is clearly narrower than the colour distribution of the GCs. 
We use the B—I broadband colour to derive a mean metal- 
licity, as it offers the largest baseline and minimizes the 
effect of random photometric errors. From the calibration 
of Couture et al. ( |1990D for globular clusters in the Milky 
Way: 



{B - I)o = 0.375 • [Fe/H] ± 2.147 



(4) 



the mean derived metallicity is [Fe/H] — — 1.01±0.05dex, 
where the uncertainty reflects only the error in the fit of 
a Gaussian function to the colour histogram. The corre- 
sponding metallicity histogram is shown in Fig. 0. The 
clearly bimodal distribution of the galactic clusters is over- 
plotted for comparison (dashed line). We cannot exclude 
that the distribution is bimodal, but then the metal-poor 
and the metal-rich peak are closer to each other than in 
the case of the galactic system. 

To search for colour gradients, we plotted (Fig. |8|) 
the B — I colour for the candidates, discarding the zone 
of the background population. Although a least-square 
fit to a function of the form B — I — A ■ 10™ ''gives 
A{B - /)/Alog(r) = 0.12 ± 0.07, the data scatter too 
much to conclude definitely that such a gradient exists. A 
similar result was obtained by dividing the sample of clus- 
ters in circular rings and calculating the mean B — I colour 
for each ring. The slope of A(B-/)/Alog(r) = 0.11±0.07 
agrees very good with the previous fit, but it is strongly 
determined by the first and the last two points, where the 
statistic is poor. From the same plot, it is apparent that 
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Fig. 6. Colour histograms of the cluster candidates. The 
bin size is 0.05 mag in each colour. A Gaussian function 
has been fit to the histograms (dashed line, see text). The 
dot-dashed line is a Gaussian with a = 0.15. The long- 
dashed line at B — I = 1.75 (lower panel) was set to divide 
the sample in red and blue clusters. 



the clusters are on average ~ 0.2 mag bluer than the un- 
derlying galaxy light. This is a property commonly found 
also in normal early- type galaxies, where the blue and 
presumably metal-poor clusters indicate the existence of 
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Fig. 7. MetaUicity histogram of the GCS of NGC 1316, 
from its B — I broadband colour. The galactic system 
(dashed line) is overplotted for comparison. The long- 
dashed line near [Fe/H] = —1.0 divides the sample in 
mctal-rich and metal-poor clusters, using the calibration 
from Couture et al ( 1990D and B — I = 1.75 (see text and 
Fig. 0.) 



a faint metal-poor stellar (halo?) popula tion, as the case 
of NGC 1380 suggests (Kissler-Patig et. al [l997| ). However, 
the non-existence of a colour gradient is consistent with 
the finding that blue and red clusters have similar surface 
density profiles (Sect. |3.3[ ). 

Differential reddening caused by the irregular dust 
structure might affect the width of the colour distribution, 
but apparently does not produce any colour gradient. 



3.2. Angular distribution 

For all objects in our sample of clusters, a transformation 
from cartesian (x, y) coordinates to polar (r, 6) was done, 
with origin in the optical centre of NGC 1316. An offset of 
50° in 9 was applied to match the PA of the galaxy quoted 
by RC3. In this way, 9 — Q represents the direction of the 
semi-major axis (sma) of the galaxy, a, and 9 = 90° the 
direction of 6, the semi-minor axis. 

To analyse the angular distribution, we rejected ob- 
jects inside a radius of 150 pixels (corresponding to 4.3 kpc 
with /X = 31.35), where the completeness is significantly 
lower (see Fig. ^ and some clusters appear over ripples 
and dust structures. Objects outside of 450 pixels were 
also rejected, as one needs equally-sized sectors to do 
this analysis, and r — 450 is roughly the radius of the 
largest circle fully covered by our frame (see dotted line in 
Fig. 0). Only candidates brighter than V = 23.8 (the 50% 
completeness level in the entire frame) were considered. 

The data were then binned in 6. Several bin sizes from 
18° to 30° were tested, and Fig. || shows the histogram 
for a bin size of 22? 5, which corresponds to dividing the 
sample into 16 sectors. The bins were taken modulo tt for 
a better statistics, that is, we assume that the distribution 
of clusters is symmetric along the semi-mayor axis of the 
galaxy. Due to the substructures present in NGC 1316, and 
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Fig. 8. The B — I colour gradient along the projected 
galactic radius. Top: the cluster candidates and a least- 
square fit (solid line). The crosses indicate the B — I colour 
of the galaxy. Bottom: the mean colour in several rings, 
with the error bars indicating the a of the mean at that 
ring. 
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Fig. 9. The angular distribution of globular clusters for 
a radius from 150 to 450 pixels (48" to 144") down to 
V = 23.8. The bin size is 22?5 and the data were taken 
mod TT. The histogram from 0° to 180° is repeated from 
180° to 360° for a better visualisation. The dashed line 
indicates the best fit of a double-cosine function. 



from the fact that it is a merger galaxy, one could expect 
some systematic differences in the azimuthal distribution 
of the clusters between both halves. Although this is not 
observed at any bin size, the small number of counts does 
not allow us to address this question clearly. 

The histogram shown in Fig. ^ demonstrates the strong 
correlation of the globular clusters with the galaxy light. 
By fitting isophotes to the galaxy, we obtained ellipticities 
ranging from 0.27 to 0.32 and position angles from 49° to 



53°, between semi-major axes of 150 to 450 pixels (the 
same used with the clusters). 

From the least-square fit to this histogram (with fixed 
period tt), we derived PA = 63° ± 9° and an ellipticity of 
0.38 ±0.06. 

There are, however, some problems that cannot be eas- 
ily resolved with our ground-based data. As indicated, 
the detection of cluster candidates near the centre of the 
galaxy is quite poor. This is, unfortunately, the most in- 
teresting region to search for young clusters which might 
be related to a merger event. 

3.3. Radial profile 

To derive the radial surface density of GCs, we divide 
our sample into elliptical annuli, as shown in Fig. |l^. The 
annuli have a width of 100 pixels along the major axis (a) 
and start from a = 50 pixels. Due to the small number 
of objects in the periphery of NGC 1316, rings beyond 
a = 900 were given a width of 300 pixels. 

The ellipticity, position angle and centre of these r ings 
were taken from the fit of the galaxy light (see Sect. 3.2) 
and fixed for all annuli. In particular, the ellipticity e, 
defined as 1 — (b/a), where a and b are the semi-major 
and semi-minor axis respectively, was set to 0.3 and the 
PA to 50°. 

We then counted the number of globular cluster can- 
didates per unit area in each ring, down to V — 23.8, 
and corrected them with the corresponding completeness 
function. 

The results are presented in Table ||. The first column 
lists the mean semi-major axis of the ring. Column 2 gives 
the raw number of clusters down to V ^ 23.8, without 
correction for completeness. Column 3 lists the corrected 
data with their errors. Column 4 lists the visible area of 
the rings, in □'. Column 5 gives the number of candidates 
per unit area. Column 6, 7 and 8 are used to compute the 
Specific Frequency (see Sect. ^). Fig. pi] shows the radial 
profile of the clusters' surface density, before and after the 
subtraction of the background counts. 

A fit of a power-law p(r) = A-r°' , where p is the surface 
density and r the projected distance along the semi-major 
axis, gives ages = -2.04±0.20 and agai -2.03±0.02 for 
the clusters and the galaxy light, respectively. We are well 
aware that a King profile may be more adequate than a 
power function, but our purpose is to compare our result 
with previous work in other GCSs in Fornax, which quote 
only power functions. 

The similarity between both slopes indicates that the 
GCS of NGC 1316 is not more extended than the galaxy 
light. Moreover, aces = —2.04 is in good agreement 
with other GCS of " norm al" early-type galaxies in Fornax 
(Kissler-Patig et al. |l997| ). 

We divided the system &t B — I = 1.75 mag into blue 
(presumably metal-poor) and red (metal-rich) clusters, 
and searched for systematic differences in the morpholog- 
ical properties between both subgroups. Again, only clus- 
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Table 5. This table gives the result of the radial profile of the cluster candidate surface density. The first column lists 
the centre of each annulus (in pixels). Then follows the raw number counts down to = 23.8, before and after the 
correction for completeness. The fourth column gives the visible area of the annulus in □'. Column 5 lists the mean 
density of GC per □'. Columns 6 and 7 give the number of clusters down to the TOM in V, before and after the 
correction for completeness. The applied geometrical corrections are listed in column 8. Finally, the number of the 
clusters in each annulus, after doubling the counts around the TOM. Note that the corrected and total number of 
clusters for the innermost annulus (sma=100) are NOT derived using the completeness correction. Instead, they were 
estimated by extrapolating the radial profile towards the centre (see text in Sect. 0). 



sma N^^s iVIf'g" area of ring GC/area TV^S'm ^Vtom geom. factor A^annui, 

[pixels] [□'] [1/D'] 



100 


20 


58.6 ± 18.1 


1.251 


46.86 ± 14.47 


20 


200.2 ± 37.5 


1. 


400 ± 75 


200 


53 


100.3 ± 17.9 


2.502 


40.09 ±7.15 


53 


100.7 ± 18.0 


1. 


201 ± 36 


300 


68 


104.4 ± 15.2 


3.751 


27.84 ± 4.05 


69 


109.4 ± 15.8 


1. 


219 ± 32 


400 


54 


65.4 ± 10.7 


5.006 


13.06 ±2.13 


56 


68.9 ±11.0 


1. 


134 ± 22 


500 


48 


58.2 ± 10.0 


6.234 


9.34 ± 1.60 


51 


62.0 ± 10.4 


0.997 


124 ±21 


600 


37 


42.9 ± 8.5 


6.159 


6.96 ± 1.22 


37 


44.5 ± 8.8 


0.821 


108 ±21 


700 


19 


22.4 ±6.2 


5.145 


4.35 ± 1.21 


19 


22.8 ± 6.3 


0.588 


78 ±21 


900 


19 


23.2 ± 6.4 


9.440 


2.45 ± 0.68 










1200 


11 


12.7 ±4.6 


7.159 


1.77 ±0.64 










1500 


4 


5.3 ±3.2 


2.880 


1.84 ± 1.11 











ters brighter than V — 23.8 were considered, leaving 150 
blue and 166 red clusters. Our results show that both sub- 
populations are equally concentrated towards the centre 
and share the same density profile (see Fig. |l2| ). 

In contrast to the radial profile, red and blue clus- 
ters systematically differ in their azimuthal distribution, 
as shown in Fig. The redder clusters show a clear cor- 
relation with the galaxy elongation, but the bluer appear 
to be spherically distributed, resembling a halo and bulge 
population. Similar results were obtained at several bin 
sizes. 



3.4. The luminosity function (LF) 

The GCLF has been widely used as a distance indicator. 
The theoretical reasoning for an universality of the TOM 
requires an initially universal mass function of GCs in 
combination with a subsequent universal dynamical evo- 
lution of the individual clusters of a GCS (Okazaki & Tosa 



1995, hereafter OT95, see e.g. Ashman & Zepf 1997 and 



Harris 2000 for reviews). While the latter is plausible, as 
far as the average evolution is concerned, the former still 
lacks a theoretical explanation. Nevertheless, the empir- 
ical evidence for an universal mass function as a power 
law diV/dM ~ M~^-^, where dN is the number of clus- 
ters in the mass range M + dM, is considerable. GCSs 
of host galaxies of very diverse nature such as the Milky 
Way, ellipticals in the Virgo cluster (Harris & Pudritz 
1994, McLaughhn & Pudritz 1996| ) or young clusters in 



merging systems like N GC 40 38/NGC 4039 (the Antenna, 
Whitmore & Schweizcr |l995|) or NGC 7252 (Miller et al. 



However, the TOM must be a function of time, being 
fainter for younger cluster systems (OT95). This is due to 
the dynamical dissolution of clusters, which affects first 
the less massive ones and then progresses towards higher 
masses. If, for example, the red population or a part of it 
was formed in the merger, one would accordingly expect 
a fainter TOM. Therefore it is reasonable to give the LF 
separately for the total system, for the red, and the blue 
subpopulation. 



The counts are given in Table |6|. As stated in Sect. p.4| , 
only objects closer to the centre than r — 300" were con- 
sidered possible members of the GCS of NGC 1316. Again, 
r is measured along the major axis of an ellipse with ellip- 
ticity 0.3. The objects with r > 300" define the back- 
ground sample. The background area is 13.720', while 
the area of the GCS is 38.010'. We then used the re- 
sults of our artificial star experiments to derive the com- 
pleteness factors in four elliptical annuli: 16" < r < 48", 
48" < r < 80", 80" < r < 176" and 176" < r < 300". 
As can be seen from Fig. |^, the completeness corrections 
differ at least in the three inner regions. The bin cen- 
ters are given in column 1. Columns 2, 4, 6 and 8 list the 
counts in each region as explained above, while columns 3, 
5, 7, 9 contain the corresponding completeness factors. 



1997| ), all show mass functions which arc well described 



by the above power law. 



The background counts, as defined in Sect. 2.4, are given 
in column 10. The total number of clusters is listed in 
column 11. This is the sum of the clusters in each re- 
gion (corrected for completeness), minus the background 
level scaled to the area of counting. The uncertainties re- 
sult from an error propagation, where also the error in 
the background was considered. A value of 0.02 has been 
adopted as the error in the completeness factors, which is 
roughly the RMS we observed between the tests. 
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Fig. 10. The elliptical annuli used in the calculation of the 
globular cluster density and specific frequency. The open 
circles indicate the position of all GC candidates brighter 
than V = 23.8. Objects outside the annulus defined by r = 
150 and r = 450 pixels (dotted line) were re ject ed for the 
analysis of the angular distribution (see Sect. |3.2| ). Objects 
outside the dashed ellipse are considered as background. 
The scale and orientation is the same as in the Fig. |l|. 



There are two additional columns, namely the number 
of blue and red clusters. They were calculated in the same 
way and will be discussed in Sect. ||. 

Bin sizes of 0.3, 0.4 and 0.5 mag were tested, and the 
results were always in good agreement with each other. 
Finally, we have chosen 0.5 mag as our bin size from the 
appearance of the fit histogram, and the absence of undu- 
lations which are present for the other cases. 

We note, however, that our derived TOM is close to 
the limit of the observations, and the last bins are strongly 
affected by the completeness correction. Nevertheless, the 
similarity of the results using different bin sizes and cen- 
ters, and the robustness of the fit against skipping the last 
bin, is encouraging. 



o 



100 



o 

(3 




100 



Fig. 11. Top: the radial profile of the density of glob- 
ular cluster candidates. The dashed line indicates the 
value adopted for the background, corresponding to 
2.0objects/n'. Bottom: the radial profile after subtrac- 
tion of the background counts. The crosses represent the 
profile of the galaxy light, arbitrarily shifted. 




2.0 

log r (arcsec) 

Fig. 12. The radial profile of the GC surface density for 
the red (triangles) and blue population (squares). No sys- 
tematic difference is seen and both are equally concen- 
trated. 

Fig. |l^ shows the three luminosity functions in B, V, 
I. For fitting the LF we chose functions, which are of 
the form: 



3V5- 



1 



(m — mo)' 
5af 



-3 



(5) 



Gaussians give no significantly different results, as ear- 
lier i nvest igations have demonstrated (e.g. Delia Valle 
et al. 1998 ). Since our photometry is not deep enough to fit 
(7t independently, we used a fixed value of at ~ 1.1, which 



is appropriate for early-type galaxies (Harris 200C quotes 
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Table 6. The counts (in bins of 0.5 mag) used in the determination of the GCLF. Given are the ^-magnitudes (or 
B, I, respectively) of the bin centers. Then follow the raw counts of the four elliptical annuli Ni (see text) together 



with the corresponding completeness factors ft. The background counts are then listed as defined in Sect. 2.4. The 
completeness factors of the fourth annulus have been used for the background as well. The total number of clusters is 
the sum of the four annuli minus the background, normalised to the same area (see text). Also given is the number of 
clusters for the blue and the red population separately, where the separating colour was B — I = 1.75. 



V 


nY 


/i 


NY 


/2 


NY 


h 


NY 


u 


J^bkg. 


"total 




^rcd 


19.5 





0.71 


1 


0.98 


1 


1.00 


2 


1.00 


0.0 


4.0 ± 2.0 


1.0 ± 1.0 


3.0 ± 1.7 


20.0 


1 


0.71 


1 


0.98 


3 


1.00 


1 


1.00 


0.0 


6.4 ±2.7 


2.4 ± 1.7 


4.0 ± 2.0 


20.5 


1 


0.67 





0.98 


4 


1.00 


2 


1.00 


0.2 


6.9 ±3.1 


3.9 ±2.6 


3.0 ± 1.7 


21.0 


2 


0.62 


6 


0.98 


3 


1.00 


3 


1.00 


1.0 


12.6 ± 5.0 


1.4 ±3.3 


11.2 ±3.8 


21.5 


2 


0.58 


3 


0.94 


12 


0.99 


7 


0.99 


2.5 


18.8 ±7.0 


7.8 ±4.6 


11.3 ± 5.2 


22.0 


2 


0.56 


5 


0.84 


14 


0.97 


11 


0.98 


4.6 


22.2 ±8.8 


14.4 ± 6.9 


8.0 ± 5.4 


22.5 


7 


0.40 


11 


0.68 


32 


0.95 


14 


0.95 


4.2 


69.8 ± 12.5 


36.6 ±9.0 


33.3 ±8.6 


23.0 


3 


0.35 


18 


0.54 


39 


0.93 


14 


0.90 


1.7 


94.2 ± 12.9 


44.3 ±8.7 


49.9 ±9.5 


23.5 


2 


0.14 


8 


0.28 


52 


0.71 


19 


0.79 


3.7 


127.2 ± 19.9 


65.5 ± 15.6 


61.7 ± 12.1 


24.0 





0.07 


2 


0.20 


32 


0.22 


8 


0.25 


3.8 


145.1 ±38.8 


48.8 ± 24.8 


96.5 ±28.1 


24.5 





0.00 





0.00 


4 


0.03 


1 


0.03 


0.5 









B 




/i 


Ni 


/2 


Nf 


h 


Ni 


k 




-'''total 


J'bluc 


N"^ 


20.5 





0.71 


1 


0.98 


2 


1.00 


2 


1.00 


0.0 


5.0 ± 2.2 


1.0 ± 1.0 


4.0 ±2.0 


21.0 


2 


0.67 


1 


0.98 


3 


1.00 


3 


1.00 


0.2 


9.5 ±3.6 


5.4 ±3.1 


4.0 ±2.0 


21.5 


1 


0.67 


3 


0.98 


4 


1.00 


1 


1.00 


0.7 


7.6 ±4.0 


0.1 ±2.7 


7.6 ±2.9 


22.0 


2 


0.58 


4 


0.96 


5 


1.00 


4 


1.00 


0.7 


14.7 ±5.0 


7.1 ±3.9 


7.6 ±3.1 


22.5 


2 


0.58 


5 


0.89 


15 


0.98 


8 


0.99 


3.3 


23.2 ± 7.9 


10.5 ± 7.0 


12.8 ±5.3 


23.0 


5 


0.42 


9 


0.74 


24 


0.96 


11 


0.96 


4.8 


46.7 ± 11.1 


32.0 ±9.7 


13.8 ±6.5 


23.5 


4 


0.35 


11 


0.64 


31 


0.94 


17 


0.92 


3.2 


70.4 ± 12.0 


35.2 ±8.6 


35.3 ±9.1 


24.0 


3 


0.26 


14 


0.32 


50 


0.80 


15 


0.85 


2.3 


127.9 ± 17.8 


57.8 ± 16.5 


70.1 ± 13.1 


24.5 


1 


0.07 


6 


0.23 


47 


0.47 


13 


0.51 


4.8 


139.8 ±27.6 


66.7 ±33.2 


73.1 ± 16.4 


25.0 





0.05 


1 


0.09 


15 


0.06 


7 


0.05 


2.0 


290.3 ± 159.0 






25.5 





0.00 





0.00 





0.00 


1 


0.02 


0.1 









7 


Nl 


/i 


Ni 


f2 


Ni 


h 


Ni 


u 




-'^t^otal 


Niiuc 


nL 


18.5 





0.71 


1 


0.98 





1.00 


2 


1.00 


0.0 


3.0 ± 1.7 


0.0 ±0.0 


3.0 ±1.7 


19.0 





0.71 


1 


0.98 


4 


1.00 


1 


1.00 


0.0 


6.0 ±2.5 


1.0 ± 1.0 


5.0 ± 2.2 


19.5 


3 


0.67 


1 


0.98 


3 


1.00 





1.00 


0.0 


8.5 ±3.3 


4.0 ± 2.3 


4.5 ± 2.3 


20.0 





0.63 


5 


0.98 


3 


1.00 


5 


1.00 


0.5 


11.7 ±4.1 


3.2 ± 2.5 


8.5 ±3.3 


20.5 


4 


0.58 


4 


0.94 


11 


0.99 


5 


0.99 


2.8 


19.8 ± 7.4 


9.5 ±4.9 


10.3 ±5.5 


21.0 


2 


0.56 


6 


0.84 


16 


0.97 


10 


0.98 


3.8 


26.7 ±8.5 


9.3 ±5.4 


17.4 ±6.6 


21.5 


4 


0.40 


9 


0.68 


35 


0.95 


17 


0.95 


5.1 


63.1 ± 12.1 


33.0 ±9.0 


30.1 ±8.0 


22.0 


5 


0.35 


18 


0.54 


37 


0.93 


12 


0.90 


1.9 


94.9 ± 13.5 


39.0 ±9.2 


55.9 ±9.8 


22.5 


1 


0.14 


9 


0.28 


46 


0.71 


15 


0.79 


2.3 


115.0 ± 17.9 


62.0 ± 13.9 


53.4 ± 10.9 


23.0 


1 


0.07 





0.20 


38 


0.22 


15 


0.25 


4.3 


199.4 ±45.4 


130.6 ±36.4 


68.8 ± 24.6 


23.5 





0.02 


1 


0.00 


3 


0.03 





0.03 


1.4 









1.4 as a typical Gaussian dispersion, while 1.29 is the ratio 
between the Gaussian dispersion and the dispersion of a 
is-function). 

Table Usts the TOMs for the three bands B, V and 
/, additionally the TOMs for the blue and the red pop- 
ulation in all filters. The uncertainties are the formal fit 
errors accounting also for the count uncertainties in the 
individual bins, as given in Table ^. 

It is apparent that the TOMs are not convincingly 
reached in the full sample. Moreover, it is clear that we 
cannot say anything definitive about a possibly different 



TOM for the red and the blue population. Deeper pho- 
tometry is required. 

The formal distance moduli, resulting from the full 
sample, are 31.58 ± 0.18 in B, 31.47 ± 0.22 in V, and 
31.19 ± 0.17 in /, adopting as the absolute TOMs M| = 
-6.89 ± 0.10, = -7.60 ± 0.0 8 and Mf = - 8.47 ± 
0.10 (Drenkhahn & Richtler |1999| , Ferrarese et al. ^000] ). 
Table 8 compiles the presently available TOMs of other 
Fornax galaxies, published by several groups. The com- 
parison shows that the B and the V TOM match the 
average values of the literature, while the / value falls at 
the faint end. However, our photometry is not sufficiently 
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Fig. 13. The aziniuthal distribution of the red (top) and 
blue (bottom) clusters. The bin size is 22? 5. As with Fig. ^, 
the data are repeated from 180° to 360° for a better visu- 
alisation. A fit of a double-sin function is overplotted for 
the red sample. 

Table 7. The turnover magnitudes for the complete sam- 
ple, and the blue and red subsamples in B, V and /. 



Filter 


TOM 


TOMbiuc 


TOM,,d 


B 


24.69 ± 0.15 


24.66 ± 0.29 


24.82 ±0.31 


V 


23.87 ±0.20 


23.99 ±0.98 


23.85 ± 0.28 


I 


22.72 ±0.14 


23.16 ±0.47 


22.49 ± 0.20 



deep, and we are not able to assign a strong independent 
weight to NGC 1316. Leaving open the problem with the 
deviating / TOM, we can only conclude that this galaxy 
has the same distance as the other Fornax galaxies. Since 
a discussion of the distance scale uncertainty is beyond the 
scope of our paper, we adopt in the following a distance 
modulus of 31.35 mag (Richtler et al. |00^. 

The contrasting finding of Gr99, who were unable to 
constrain the TOM in spite of their deeper data, may indi- 
cate that the central region, where Gr99 performed their 
analysis, hosts a different population of objects (see the 
discussion). 



4. Specific frequency 

The specific frequency of a GCS is defined as the num- 
ber of clusters per unit luminosity of the parent galax y 
referred to My = — 15 (Harris & van den Bergh 1981 ): 
Sn = N ■ ioO-4 (Mv+i5)^ .^here N is the total number of 
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Fig. 14. The luminosity functions for B, V and /. The 
bin size is 0.5 mag. The last bin (open circle) was not used 
in the fit, as its center is beyond the limiting magnitude. 
The smooth curve is the fit of a is function with cr = 1.1. 



clusters and My the absolute ^-magnitude of the host 
galaxy. 

The basis for calculating N is Table |5[ We counted the 
number of globular clusters in each elliptical ring down to 
the TOM in V (see Sects. U and 

The data in each ring were then corrected by the corre- 
sponding completeness function. When the rings were not 
fully covered by our frames (see Fig. |l^), a geometrical 
correction factor was used to scale the counts to the total 
area. The results are listed in Table |. The raw number of 
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Table 8. Published TOMs for galaxies in the Fornax clusters, using the GCLF. References are as following: 1. Kohle 
1996t 2. Blakeslee & To nry|1996| ; 3. Elson et al. |1998| 4. Delia Valle et al. 



et al 
et al 



1999; 7. Richtler et al. 



1992 



199g; 5. Ostrov et al. 1998: 6. Grillmair 



Galaxy 



TOM(B) TOM(V) 



TOM(/) 



Reference 



NGC 1344 








23.80 


± 


0.25 








2 


NGC 1379 


24.95 


± 


0.30 


23.65 


± 


0.26 


22.77 


± 


0.44 


1,3 


NGC 1380 


24.38 


± 


0.09 


24.05 


± 


0.25 








2,4 










23.67 


± 


0.11 








4 


NGC 1399 


24.59 


± 


0.10 


23.83 




0.15 








2,6 










23.90 


± 


0.08 


22.27 


± 


0.05 


1 










23.71 


± 


0.12 








5 


NGC 1404 


24.86 


± 


0.21 


24.10 


± 


0.20 








6,7 










23.92 


± 


0.20 








2 


NGC 1374 








23.44 


± 


0.13 


22.49 


± 


0.12 


1 


NGC 1427 








23.63 


± 


0.19 


22.22 


± 


0.13 


1 



100 




Fig. 15. The GCLF of the blue and the red populations, 
in the V band. 



counts are listed in column 6. Column 7 lists the counts 
after the completeness correction. The geometrical factor 
is given in column 8. To calculate the total number of clus- 
ters in each annulus (column 9), we assumed the luminos- 
ity function to be symmetric around the TOM and dou- 
bled our counts. The total population of clusters (which 
includes background objects) was derived as the sum of 
the counts listed in column 9. The last three annuli were 



not included, as the counts cannot be distinguished from 
the background. Finally, to estimate the number of clus- 
ters in the central region (first annulus and inner ellipse) 
we extrapolated the data to follow the slope of the radial 
profile defined by the outer bins. This is because of the 
poor detection of clusters in this region, which would make 
the completeness correction too uncertain (see Fig. ||). 

Several authors argue that tidal disruption and dy- 
namical effects are more effecti ve near the centre (for a 
discussion see e.g. Gnedin 1997). Our purpose is to show 
that even if the central GCs were not affected by these de- 
struction mechanisms, the specific frequency is still very 
low. 

We then arrived at a mean density of 160 ± 
30objects/n' down to the TOM for the innermost bin and 
inner ellipse. Thus, for this region of area 1.407n', 225±42 
clusters brighter than the TOM are expected. As with the 
outer bins, we doubled our counts around the TOM to 
obtain iVa„nuius = 450 ± 84. 

Summing up the results for each annulus gives 1314 ± 
240 (see column 9 of Table ^l]). From this, one only 
needs to subtract the background counts, estimated to be 
2.0 objects/n'. As we considered the geometrical incom- 
pleteness of the last three bins (see Fig. |l^ and column 8 in 
Table |ll|), the total area of counting is 35. 1860' including 
the nucleus. Therefore, about 141 background objects are 
to be subtracted. Again, the counts were doubled around 
the TOM. Finally, we estimate Nt = 1173 ± 240 as the 
total population of clusters in NGC 1316. 

With n 
magnitude my 



31.35 (Richtler et al. 2000) and an apparent 
5.53 (de Vaucouleurs et al. 1991), we 



obtain Sn = 0.9 ±0.2, where the uncertainty results from 
error propagation and has the form: 



a = ioO-4-(Mv+i5) ^^2^ ^ (0.4 • iV • In 10)^ 



'Mv 



(6) 



This is an extremely low value for a bright early-type 
galaxy. As we discuss below, it may result from the fact 
that the dominant light contribution to NGC 1316 comes 
from an intermediate age population instead from an old 
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one. A similar interpretation has already been made by 
Hilker & Kissler-Patig ( |1996D in the case of NGC 5018. 

Furthermore, if one assumes that, due to dynamical 
effects, the surface density of the cluster population does 
not follow the power law, but tends to remain approxi- 
mately constant in the two inner annuli, then the number 
of clusters drops to about 150, and the Sn is decreased to 
0.7. 

Nevertheless, we cannot exclude the possibility of mea- 
suring an erroneous background density. As can be seen 
from the radial profile (Fig. |ll|), the background level is 
reached at r ~ 300 arcsec, where the average density is 
2.0objects/n'. It is interesting to investigate what are the 
effects of varying this number by la, as it enters directly 
into the Specific Frequency and hence, into the estimated 
age of the merger. This is discussed in Sect. |5.2|. 



5. Discussion 

5.1. Subpopulations of clusters 

The detection of intermediate-age clusters by GoOO clearly 
shows that the cluster system is inhomogeneous in terms 
of age. Also the strikingly different angular distributions of 
the blue and red clusters argues in favour of the existence 
of subpopulations. One would naturally assume that the 
red and probably metal-rich distribution, or part of it, has 
been formed in the merger. A popular scenario attributes 
the observed bimodal colour distributions in many early- 
type galaxies as being due to cluster formation in enriched 
gas during a merger (Ashman & Zepf 1992, 1997). This 



bimodality is not seen in NGC 1316. However, we think 
that this is not of great significance, given the following 
points. 

A bimodal metallicity distribution as a result of clus- 
ter formation during a merger probably needs special con- 
ditions. If, for example, the merger components were two 
spiral galaxies with a certain range of metallicities present 
in the gas content of their disks, then one naturally expects 
cluster formation with a similar range of metallicities, per- 
haps weighted to higher metallicities due to chemical en- 
richment during the merger. Any bimodal metallicity dis- 
tribution, which has been present before the merger (and 
which does not stand in relation with the merger, like the 
bulge and disk clusters in the Milky Way) , could well have 
been blurred. 

Moreover, a metal-rich intermediate-age population of 
clusters, formed in the merger, for example with an age of 
2 Gyr, would have B — V colours which are about 0.2 mag 
bluer than their old, metal-rich counterparts. According to 
the above relation between B — I colour and [Fe/H], this 
corresponds to a metallicity difference of 0.5 dex. The peak 
at —0.8 dex in Fig. could well be explained in this way, 
comprised perhaps by clusters with solar metallicity and 
those found by GoOO. 



5.2. The Specific Frequency 

Our result of 0.9 ± 0.2 for the specific frequency of the 
GCLF of NGC 1316 confirms with a larger field what has 
been found by Gr99 using HST-data. This is an extremely 
low value for a bright, early- type galaxy, even more aston- 
ishing as one would expect to have many clusters formed 
during the merger process. 

One possible explanation for this low value would be 
to assume that the merger caused, whether or not star 
clusters have been formed, a period of star formation, 
which nowadays constitute an intermediate-age popula- 
tion, which is brighter than an old population without 
strikingly differing in integral broad-band colours. The 
contribution of such a population could be visible spec- 
troscopically by stronger hydrogen lines than normal. 

We assume that a "normal" specific frequency of early- 
type galaxies is 4. The question then is: How young must 
NGC 1316 be to lower that value to our value of 0.9? 



Worthey (1994) offers a tool for synthesizing stellar pop- 
ulations. Assuming a metallicity of -1-0.5 dex, we get from 
his "model interpolation engine" an age of 2 Gyr to ac- 
count for the magnitude difference of 1.6 mag in compari- 
son to an old population of 15 Gyr. 

The dominant error in this result is the uncer- 
tainty in the background density. If it were, instead of 
2.0objects/n', 4.0objects/n', the total number of clusters 
decreases to Nt = 1033 and Sn = 0.8. On the other 
hand, a density of O.Oobjects/D' implies Nt — 1314 and 
Sn — 1.0. Thus, a total error of ~ 0.3 results for Sn- This 
means, from Worthey's models, an error in the age of the 
merger of ~0.8Gyr. 

Our value of 2 Gyr is in good agreement with the 
luminosity-weighted mean stellar age of NGC 1316 given 
by Kuntschner & Davies (1998), who measured line 
strengths for early-type galaxies in the Fornax cluster. It 
also matches the merger age of NGC 1316 estimated by 
Mackie & Fabbiano ( |1998D and the globular cluster ages 
of GoOO. In this context, one also should note that the 
occurence of two SNe la is a further indication for the ex- 
istence of an intermediate-age population, to which many 



SNe la seem to belong (e.g. McMillan & CiarduUo |996|). 

The question now arises whether specific frequencies of 
GCS in general could be influenced, if not dominated, by 
the existence of intermediate-age populations. One notes 
that among t he sa mple of early-type galaxies, for which 
Trager et al. (2000) measured age-sensitive line strengths, 
a few striking examples like NGC 221 or NGC 720 exist, 
which have a young component and whose specific fre- 
quencies are quite low for early- type galaxies. Another ex- 
ample is the merger candidate NGC 5018, for which Hilker 
& Kissler-Patig ( 1996) found a low specific frequency, sim- 
ilar to that of NGC 1316, which they also explained by the 
presence of a young population. 

A strong argument against the conjecture that this ef- 
fect is dominating, is the fact that the specific frequency 
among early-type galaxies has a positive corre lation with 
the environmental galaxy density (West |l993 ). If an en- 
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hanced number of mergers in a high density environment 
caused the existence of intermediate-age populations in 
many early- type galaxies, one would expect a negative cor- 
relation, if the brightness of the host galaxy significantly 
influenced the specific frequency. However, if the efficiency 
of globular cluster formation increases with the star for- 
mation rate, as the work of Larsen & Richtler (2000) in- 
dicates, then a positive correlation, as observed, is under- 
standable. 

Larsen & Richtler ( 200C ) performed a survey of young 



clusters in a sample of 21 nearby face-on spiral galaxies. 
They found that the "specific luminosity" Tl{U), which 
measures the ratio of the summed-up U-luminosity of clus- 
ter candidates to the total U-luminosity of the host galaxy, 
correlates well with indices measuring the star formation 
rate. Apparently, a high star formation rate causes also a 
higher cluster formation efficiency. 

The final specific frequency after a merger therefore 
depends on many circumstances, among them the star 
formation rate, the amount of gas available for star for- 
mation, the specific frequency of the merger components, 
etc. 

Can we expect that the merger enhances the specific 
frequency of the outcoming configuration? GoOO describe 
the evidence that makes NGC 1316 a candidate for an ad- 
vanced disk-disk merger system. So let us assume that two 
spiral galaxies of equal brightness My have merged. 

If we set for the visual luminosity L ~ iq-oa-Mv then 
the following relation holds: 



iVm = 2 • 10"'' ■ Sn-L 



(7) 



with iVtot the total number of clusters, N^^ the number 
of cluster formed in the merger, and Sn the specific fre- 
quency of the merging spiral galaxies. Further: 



Ltnt — 2 ■ L, 



dim 



in 



(8) 



where Ltot is the actual luminosity of NGC 1316, idim the 
actual luminosity of the former spirals, dimmed for about 
3 Gyr, after star formation ceased, and is the lumi- 
nosity of the stellar population formed in the merger. A 
typical value for Sn is 0.5 (Harris 2000| ). N^^ is not accu- 
rately known, but a rough estimate comes from the fact 
that GoOO found 3 clusters of masses comparable with 
LO Cen, so a lower boundary is 300 clusters. With these 
numbers, the luminosity of the hypothetical spirals, at 
the time of merging, turns out as 9 • 10^, corresponding 
to My — —22.4. The dim factor during the next 3 Gyr 
is about 2, based on a simulation with Worthey's "model 
interpolation engine". If Ltot is 1.32 • 10^, then L,„ turns 
out as 4.2 • 10^, corresponding to My = —21.6. The stel- 
lar population formed in the merger will dim for 1.5 mag 
during the next 12 Gyr. The further fading of the spiral 
will be roughly by factor of 5, corresponding to the ratio 
of 10 between the M/L values of a spiral and an elliptical 
galaxy. Thus the total luminosity after aging will be about 
3-10^, which means My — —21.2 and a specific frequency 
of 4. 



The numbers used in this exercise are uncertain, but it 
shows that the higher specific frequency of elliptical galax- 
ies with respect to spirals can indeed plausibly be under- 
stood as the outcome from disk-disk mergers. 

The example also shows that a merging scenario of an 
elliptical galaxy with a spiral leads to such high merger 
luminosities, that the gas mass consumed by the star for- 
mation in the merger cannot reasonably be brought in by a 
spiral galaxy. For instance, an elliptical with My = —21.5 
and specific frequency 3 corresponds through the above re- 
lations to a spiral galaxy with My = —21.2 and a merger 
brightness of My = —22.4. A population of this bright- 
ness with an age of 3 Gyr has a mass of almost 8 • 10^° Mq, 
which cannot be provided by any ordinary-sized spiral. 

It is, however, doubtful, whether high specific frequen- 
cies (6 and more) can be created by the merging of two 
or more spiral galaxies (Harris |200C ). It is more likely 
that dwarf galaxies must be involved if cluster systems 
like t hose o f NGC 1399 and M87 are to be formed (Hilker 
et al. 



199C) 



5.3. The Luminosity Function 

The importance of the LF of NGC 1316 lies in the fact that 
it could provide an accurate distance to the two SNe la 
SN 1980N and SN 1981D and thus help to cahbrate the ab- 
solute SN la luminosity, as has b een t he case for SN 1992A 
in NGC 1380 (Delia Valle et al. |l998| ). 

However, due to its presumed merger history, one 
might suspect that the CCS of NGC 1316 is contami- 
nated by younger clusters, which very probably should 
be searched for among the red cluster population. Broad 
band colours are not very sensitive: a cluster becomes red- 
der in y — / by only 0.2 mag, evolving from an age of 2 Gyr 
to an age of 15 Gyr. 

Such a contamination would also affect the turnover 
magnitude. OT95 investigated the evolution of a GCS 
considering the dynamical and photometric evolution of 
individual clusters. According to their results, the TOM 
would brighten by about 0.6-1.0 mag, when a GCS evolves 
from an age of 2 Gyr to an age of 15 Gyr. Therefore one ex- 
pects the TOM of the red cluster population to be fainter 
than that of the blue population, if the red clusters are 
indeed younger. 

This has to be proved/disproved with deeper data than 
those at our disposal. Since the TOM is not reached in 
either among the red population or the blue population, 
its numerical value returned by the fit is rather uncertain 
and hence we cannot make any secure statement. 

From the models of OT95, we would expect to see the 
TOM approximately at V = 24.5 or B = 25.5 in case that 
the majority of the red clusters has been formed in the 
merger. 

Gr99, in a small field near the centre of NGC 1316, 
obtained _B-photometry with HST to faint magnitudes. 
Their results are not easy to interpret. Their actual 
background-corrected and completeness-corrected counts 
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do not rule out a turnover at B ~ 25.5, but no turnover is 
indicated in their /-counts at / = 23.7, where one expects 
it to be due to the mean B ~ I colour of 1.8. If they cor- 
rect for the contribution of "dust-obscured" objects, the 
resulting LF is completely dominated by the unseen pop- 
ulation and there is a steady increase in both B and / 
down to the faintest magnitudes. 

Gr99 conjecture that they see a population of open 
clusters, which is plausible if the star formation near the 
center has continued for some time after the merger. One 
can imagine that most of the globular clusters formed dur- 
ing the merger, when the star formation rate was high, but 
less massive open clusters continued to form at lower rate, 
which would blur the luminosity function of the globular 
clusters. 

No detailed picture follows from the present data, but 
it is not excluded that a large fraction of the red globular 
clusters formed during the merger. Another (unknown) 
part may have been brought in by the merger components. 
If two spiral galaxies have merged they could have brought 
in their metal-rich bulge populations. The blue globular 
clusters would then be the common halo globular clusters 
of the two merger components. Star formation near the 
centre continued and younger star clusters were added to 
the cluster population formed during the merger. In that 
case, we would expect the outer population of red clusters 
to exhibit a TOM in the range 0.6 to 1 mag fainter than the 
blue cluster population, which could be shown by deeper 
photometry than we have. 

In this scenario one should only regard the LF of the 
blue cluster population as appropriate for distance esti- 
mation. 



5.4. A remark on the supernovae in NGC 1316 

Two SNe of type la, SN 1980N and SN 1981D, appeared in 



NGC 1316 (Hamuy et al. |l99lD . The distance determined 
by its GCLF can thus contribute to the calibration of the 
absolute peak brightness of SN la. There are only a hand- 
ful of nearby early- type host galaxies for la supernovae, 
so NGC 1316 occupies an important role. 

However, our data are not deep enough to measure the 
distance via the GCLF so accurately that SN 1980N could 
be independently compared with SN 1992A in NGC 1380, 
the third SN la which appeared in the Fornax cluster, and 
which has a very well measured GCLF (Delia Valle et al. 



1998| ). SN 1981D has a hght curve of lower quahty. 

As Tables g and | show, the TOM of the GCS of 
NGC 1316 fits at least well to the other TOMs of early- 
type Fornax galaxies. This may be taken as an argument 
that NGC 1316, in spite of its outlying location, is in the 
same distance as the core of the Fornax cluster. 

Indeed, after the appropriate corrections for decline- 
rate and colour (e.g. Hamuy et al. 199(:, Richtler et al. 
2000| ), the apparent peak brightness of SN 1980N agrees 
excellently with that of SN 1992A, which is a much 
stronger argument for their equal distance. 
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